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Abstract 

A realistic computer model for simulating isothermal and non-isothermal phase transformations proceeding by homoge- 
neous and heterogeneous nucleation and interface-limited growth is presented. A new treatment for particle size effects on 
the crystallization kinetics is developed and is incorporated into the numerical model. Time-dependent nucleation rates, 
size-dependent growth rates and surface crystallization are also included. Model predictions are compared with experimental 
measurements of DSC/DTA peak parameters for the crystallization of lithium disilicate glass as a function of particle size, 
Pt doping levels, and water content. The quantitative agreement that is demonstrated indicates that the numerical model can 
be used to extract key kinetic data from easily obtained calorimetric data. The model can also be used to probe nucleation 
and growth behavior in regimes that are otherwise inaccessible. Based on a fit to data, an earlier prediction that the 
time-dependent nucleation rate in a DSC/DTA scan can rise above the steady-state value at a temperature higher than the 
peak in the steady-state rate is demonstrated. 


1. Introduction 

As was demonstrated in a companion paper in this 
volume [1], * differential scanning calorimetry (DSC) 
and differential thermal analysis (DTA) provide con- 
venient and rapid methods for studying phase transi- 
tions under isothermal and non-isothermal condi- 
tions. By carrying out properly designed experi- 
ments, the relative importance of surface and volume 
nucleation can be assessed and some indication of 
the temperature range of significant nucleation can 
be obtained. These calorimetric techniques suffer, 
however, from the paucity of methods for quantita- 


* Corresponding author. Tel.: -I- 1-314 935 6228; fax: -I- 1-314 
935 6219; e-mail: kfk@wuphys. wustl.edu. 


tive data analysis. Analysis techniques are typically 
based on Kissinger-type methods [2,3], which as- 
sume an Arrhenius temperature dependence for the 
transformation kinetics. Since the kinetics are dic- 
tated partially by the nucleation rate, which is defi- 
nitely non-Arrhenius, the use of these methods is 
questionable when applied to most first-order phase 
transformations. Experimental [4,5] and theoretical 
[6] investigations have shown that such techniques 
can be used with some confidence in transformations 
involving growth only. There the activation energy 
obtained is approximately equal to the activation 
energy for growth. Since this requires that either the 
primary nucleation be heterogenous and that the sites 
are saturated, or that the peaks in the nucleation and 
growth rates be widely separated, these methods 
have limited applicability. Also, since growth is 
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strictly non-Arrhenius, the activation energy obtained 
will depend on the temperature range over which the 
transformation took place. This can depend critically 
on the number of quenched-in nuclei or the rate at 
which the scan takes place through the nucleation 
zone, where there will be no accompanying DSC or 
DTA signal. Clearly, better methods of analysis that 
can take such effects into consideration are required. 

We have developed a realistic numerical com- 
puter model for simulating isothermal and non-iso- 
thermal phase transformations proceeding by homo- 
geneous and heterogeneous nucleation and 
interface-limited growth. Unlike earlier models, the 
nucleation rate need not be steady-state; transient 
nucleation which is dominant in the formation and 
crystallization of many silicate and metallic glasses, 
is included naturally. Further, no ad-hoc assumptions 
are made of what constitutes a growing cluster con- 
trasted with a nucleation embryo. Cluster evolution 
is followed directly from the smallest embryo to a 
macroscopic cluster. Previous comparisons between 
the model predictions and experimental non-isother- 
mal DSC and DTA data for a polymorphically devit- 
rifying glass, lithium disilicate, are in good agree- 
ment. The computer model has been used to demon- 
strate that all existing methods for analyzing non-iso- 
thermal DSC data are flawed [6] and to evaluate a 
recently proposed method for estimating nucleation 
rates from DSC and DTA data [7], 

Those calculations were made assuming volume 
nucleation and growth at constant composition in 
large samples. As was discussed in our companion 
paper, however, an ability to model finite-size effects 
and surface nucleation and growth can often be 
critical for a quantitative analysis of experimental 
data. Here we extend the numerical model to include 
these features. It should be emphasized that this is 
not the first attempt to include finite-size effects and 
surface crystallization in the analysis of phase trans- 
formation kinetics. Surface crystallization has been 
modeled for both metallic [8] and silicate [9] glasses. 
These earlier treatments are followed for modeling 
the surface crystallization here. Finite-size effects on 
isothermal transformation kinetics have also been 
considered previously by Weinberg [10, 1 1]. The ex- 
pressions developed there, however, are mathemati- 
cally complex, making them difficult to extend to 
non-isothermal transformations. Our simpler expres- 


sion is more easily used to analyze isothermal trans- 
formation kinetic data and is readily adapted to more 
complicated situations, such as the multi-step 
isothermal and non-isothermal annealing treatments 
considered here. The validity of our expression has 
been established by a quantitative comparison with 
computer simulations of phase transformations on a 
lattice of discrete, non-interacting, spins. 

To test the extended computer model, we com- 
pare the theoretical predictions of DSC/DTA data 
taken under conditions of non-isothermal heating 
with the experimental studies of the devitrification of 
lithium disilicate glass discussed in our companion 
paper. An important result of this modeling is the 
confirmation of the importance of time-dependent 
nucleation in determining DSC/DTA behavior. In 
particular, the agreement between model predictions 
and experimental data for isothermal and non-iso- 
thermal anneals in the temperature range of signifi- 
cant nucleation supports our previous model predic- 
tion that the time-dependent nucleation rate in a 
DSC/DTA heating scan peaks at a temperature 
higher than the steady-state rate and rises signifi- 
cantly above the steady-state value there. 


2. Description of the numerical model 

The basic numerical model and the new exten- 
sions to include particle size and shape are discussed 
in this section. Since the basic model has been 
discussed extensively elsewhere [6,7,12-14], only 
the most salient points are reviewed here. 

2./. The basic model for an infinite sample 

The classical theory of nucleation [12] is assumed, 
taking direct account of time-dependent nucleation 
rates. Assuming spherical clusters, negligible stress 
effects and a sharp interface between the cluster and 
the parent phase, the reversible work of formation 
for a cluster of n molecules can be written as 

W n = n\G f T (367T ) !/ ? 'Tr / ^n 1/ \r . ( 1 ) 

Here, AG' is the Gibbs free energy per molecule of 
the new phase less that of the initial phase, v is the 
molecular volume, and a is the interfacial energy 
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per unit area. Clusters are assumed to evolve slowly 
in size by a series of bimolecular reactions, leading 
to a time-dependent cluster distribution described by 
a system of coupled differential equations of the 
form 

d N nt 

= +*.a + ] 

+ a Ui.A+p ( 2 ) 


duration 8/ = 57/0, where 0 is the scan rate in 
degrees per second and 57 is the temperature step 
size, allowing the nucleation rate to evolve in each 
interval as described by Eqs. (2) and (4). At the end 
of the interval, the sizes of the nuclei generated in 
previous intervals are calculated using the average 
growth rate, w(r), given by Eq. (5). The extended 
volume fraction transformed, .v c , is calculated at the 
end of each interval: 


where A/, is the number of clusters of size n at time 
t. The rate of monomer addition and loss, &,} and 
k~ respectively, are obtained from the diffusion 
coefficient in the parent phase and the atomic jump 
distance (see Refs. [12,13] for details). The nucle- 
ation rate at a cluster size w, I n r is the time-depen- 
dent flux of clusters past that size and is given by 

= -*„ + (3) 

The coupled differential equations are solved us- 
ing a finite difference method (see Kelton et al. [13]) 
in which the time is divided into a large number of 
small intervals, 8/, and the number of clusters of size 
n at the end of the interval, N nt + hr is calculated 
from 


d N 


N n.t+St ~ N nj + 


d t 


(4) 


where dN n Jdt is given in Eq. (2). 

The behavior of clusters up to some large size is 
calculated directly by these methods. The growth of 
clusters larger than this size is computed using an 
approximate expression for the growth rate as a 
function of cluster radius, w(r), that was derived 
previously [14]: 


16D 

«('•) = c— 


3D 

4 ir 


i/3 

sinh 


2 k B T 



(5) 


Here D is the diffusion coefficient in the parent 
phase, A is the jump distance, A G v is the Gibbs free 
energy per unit volume, 7 is the temperature, k B is 
Boltzmann’s constant, and C is a constant that is 
dependent on the mode of growth [15]. 

To simulate glass formation and DSC non-isother- 
mal scans, the rate of volume fraction transformed 
under non- isothermal conditions is computed by di- 
viding the time into a series of isothermal scans of 



( 6 ) 


where is the number of nuclei generated in the 
interval /, r i } is the time-dependent radius of those 
nuclei, and V 0 is the sample volume. Assuming that 
the nucleation and growth occur randomly in space, 
the actual volume fraction, v(r), is computed assum- 
ing Johnson-Mehl-Avrami-Kolmogorov (JMAK) 
statistics [16-18]: 


x(t) = 1 — exp( —x e ). 


(7) 


Assuming a linear relation between the rate of change 
of the enthalpy, dH/dt , and the rate of volume 
fraction transformed, the DSC trace is computed as a 
function of temperature or time as 


DSC signal a 


x(r, + g 7 )-A( 7;.) 

St 


( 8 ) 


2.2. Model extensions 

The model discussed in the previous section as- 
sumes (1) that the sample is large compared with the 
size of an individual grain of the daughter phase, (2) 
that nucleation occurs randomly in space and time 
and (3) that the growth is three dimensional with the 
grains of the daughter phase having a low aspect 
ratio. Devitrification in real samples, however, often 
occurs under quite different conditions. Surface crys- 
tallization, for example, is common in silicate glasses, 
although internal, or volume, nucleation and growth 
usually dominates in bulk glasses since the surface 
constitutes a small fraction of the sample. Samples 
used for calorimetric studies, however, are generally 
fine-grained powders (typically 50 to 400 (am in 
diameter) with a very large surface to volume ratio. 
Surface growth, if it occurs, will then contribute 
more substantially to the transformation signature. 
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Also, since a cluster nucleates within a particle, all 
space is not available for growth; the region outside 
the particle is excluded, thus violating a fundamental 
assumption of the JMAK statistics. In this section, 
we present approximate methods for dealing with 
these effects in isothermal and non-isothermal trans- 
formations. Their successes and limitations are dis- 
cussed in Section 3. 


2.2./. Finite size (FS) effects 

The probable number of nuclei appearing in a 
smalt particle of volume F p , given the number ex- 
pected in a mole of volume V mok , is N p ~ 
jY nui | 0 ( Vp/ V nu) | c ). Expressed in terms of the nucleation 
rate, / mole , then the average number of nuclei per 
particle appearing in a given time interval, A r, is 


N P = A/. 

V. 


( 9 ) 


The distribution of nuclei among the particles of 
the ensemble, given N p average nuclei per particle, 
is calculated from the Poisson distribution: 


N"r 

/> (« P ) =e ~ ,p 7 ’ ( l0 ) 

"p ! 

where // is the number of nuclei per particle. To 
minimize the calculation time in the computer mod- 
eling, slightly different approaches were taken de- 
pending on the size of I V . For N p < 300, Eq. (10) 
was applied for each value of n p between 0 and 5 A p , 
a value sufficiently large that P(5N ? )< 1. For 300 
< 10000, the particles were partitioned into 
bins of width A (set by the constraint that the total 
number of bins was 1500) and Eq. (10) was used to 
compute Pin ), taken as the average value over the 
bin (i.e., (Pin))). For N p > 10000 the distribution 
was sufficiently narrow that statistical fluctuations 
among particles could be ignored. 

Growth presents different problems. Given a sin- 
gle nucleation event within a particle of radius R , a 
cluster will grow' w r ith a cluster-dependent growth 
velocity given by Eq. (5). If the cluster were to 
nucleate at the exact center of the particle and if 
there were no other growing clusters wuthin the 
particle, the transformation would proceed unim- 
peded until the particle was completely transformed. 
For clusters nucleated away from the center, how- 




Fig. 1. Schematic Figure illustrating the case of a spherical nucleus 
growing in a finite-sized particle. 


ever, the transforming region will intersect the parti- 
cle surface prior to the complete transformation of 
the particle. Now' the assumption of unimpeded 
growth requires that some of the transformed volume 
fall outside the boundary of the particle (labeled as 
V (HJ , in Fig. 1). Clearly this region cannot contribute 
to the actual volume fraction transformed. To model 
such finite size effects, we assume that (1) the 
clusters form randomly in space and time and grow 
independently of each other, (2) the number of nuclei 
produced in a particle is proportional to the volume 
of the particle (Eqs. (9) and (10)) and (3) growth 
outside the particle is prohibited. Since finite com- 
puter resources make it impossible to calculate the 
positions of all clusters in all particles and to follow 
their growth directly, a statistical description is de- 
veloped which gives the average growth volume 
outside a particle in terms of the cluster and particle 
radii, r c and R respectively. This volume, V^ ul (r r , 
/?), is then subtracted from the extended volume 
transformed before applying the JMAK statistics (Eq. 
( 7 )). 
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For simplicity, we assume spherical particles of 
radius R . Given a cluster of radius r c with a center 
at some position £< R, the volume of the cluster 
falling outside the boundary of the particle is 

V'( /?, r c ) = — r£( 1 - cos 0 2 )~(2 + cos 0 2 ) 

7T 

/?■(! — cos 0, )~(2 + cos 0, ) , 


R~-r~-C 


C 2 d£ 


[K 



7T 

J H-r c 

1 - 

l 2R C J 


R--r{ + r\ , 

X 2+ dc 

2 RC * * 


where 0] and 0 2 are defined as shown in Fig. 1. 
Three cases are relevant to the calculations pre- 
sented, corresponding to different positions and rela- 
tive sizes of the growing cluster and particle: (i) 0, 
and 0 2 < 90°; (ii) 0, < 90° and 0, > 90°; and (iii) 
0, and 0 2 > 90°. An application of the law of cosines 
demonstrates that for all three cases V'(R , r c ) is 
given by Eq. (11). These results can be extended 
directly to the case of clusters growing in a large 
number of particles of equal size, constituting a 
powder sample, such as is often used in DSC or 
DTA experiments. For a cluster w ith a radius smaller 
than the particle in which it is growing to have 
predicted growth outside the particle, the cluster's 
center must be located within a shell of inner radius 
R - r c and outer radius R. Assuming the internal 
nucleation to be equally probable in all regions, the 
probability that the center of a given cluster lies 
between a radius £ and f -Fd£ in the particle is 

4tt£ 2 df 2 

= 7 = — rd£. (12) 

v ’ (4i r/?-’)/3 R- v ’ 

The average value for the volume lying outside the 
particle, V om ( R* r c ) can be obtained by integrating 
V r { /?, r c ) over all possible locations of the center of 
the growing cluster in the particle, weighing by the 
factor in Eq. (12). Since clusters with centers farther 
than r c from the particle surface grow entirely within 
the particle, they will not contribute to the volume 
outside the particle: hence the cluster center loca- 
tions of interest lie within the bounds R — t\- to /?, 
givine 


KJ R < r ( ) 


Because the weighing increases strongly with in- 
creasing r, the lower limit can be extended to r = 0 
without significant error, leading to 

7T/v 4 18 R 2 — rr 

KA R ->\) = — rr— • (14) 

24 a 

This formulation is correct for all cases where 0 < r c 
<2 R. Clearly, for those cases where the growing 
cluster has a radius larger than that of the particle 
(R < r c < 2/?), the lower limit of the integral must 
be set to zero to avoid double counting. Values of 
V oul (R, r c ) computed from Eq. (14) agreed quantita- 
tively with the results from a direct computer calcu- 
lation of the average cluster volume lying outside the 
particle, for a set of clusters of given size that were 
located randomly in an ensemble of 1 0000 spheres. 
For more than one growing cluster, the average value 
for V ini[ (R , r c ) per particle must be obtained by 
summing the separate contributions of all growing 
clusters: 

Kl,= E 'V). (1?) 

all clusters 

where V oul (/?, r c ) is the average volume lying out- 
side the particle for a cluster of radius / t . F° r the 
purposes of the calculation, then, the volume trans- 
formed, taking into account that which lies outside 
the particle but not accounting for overlap of clusters 
within the particle, is 

E f^c-vu *.'•<-))■ ( 16 ) 

all clusters ' “ t 


r c 

h r,- 


R - - - £ 


and the ‘corrected’ extended volume fraction trans- 
formed. a',, is 

V V 

Vl = 7T = (17) 
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A powder will actually contain a distribution of 
particle sizes, which can be incorporated easily into 
the calculation. A series of equations having the 
same form as Eqs. (16) and (17) is used, one for 
each particle size, /?,, giving 

V'= L ( l8 > 

all clusters, \ ^ ' 

all particles 

The corrected extended volume fraction transformed 
is changed in a similar manner: 

V' 

.v'. = t . (19) 

L E B . l|W ucWs(4ir« 3 )/3 

The computer simulation of the transformation now' 
proceeds as for the infinite-sized particles, using Eqs. 
(16) or (18) to compute the extended volume fraction 
transformed in each iterative step and computing the 
corrected extended volume fraction transformed us- 
ing Eqs. (17) or (19). Overlaps between different 
clusters are assumed to obey the standard JMAK 
statistics (Eq. (7)). giving the volume fraction trans- 
formed as 

.v — 1 ~ exp( — a',) . (20) 

The DSC/DTA signal is then computed using Eq. 

(8). 

It is important to emphasize that this method for 
incorporating finite particle sizes still uses the JMAK 
method to account for cluster overlap. That analysis 
hinges on the assumption that the nudeation and 
growth occurs randomly w’ithin the sample and that 
the size of each transformed cluster is small com- 
pared with the sample volume. Because of this last 
point, the analysis developed in this section should 
be valid when the number of nuclei in each cluster is 
large. A high nudeation rate during the scan or a 
large number of quenched-in nuclei is therefore re- 
quired. Unfortunately, this condition is frequently 
not true in silicate glasses, given the relatively low 
nudeation rates. A further extension of this analysis 
is discussed in the next section. 

2.2.2. Corrections to finite-size calculation for small 
numbers of nuclei 

For small particles, when the number of nuclei is 
less than 5 to 10 per particle, non-isothermal calcula- 


tions made following the approach discussed in Sec- 
tion 2.2.1 predict two DSC peaks, which is in con- 
flict with the experimental data, showing only one 
peak for all particle sizes. This is an artifact arising 
from a failure of Eqs. ( 1 3)— (20) to properly account 
for the excluded volume when the volume of each 
transforming region becomes a significant fraction of 
the particle volume. This was confirmed recently 
from a comparison between the predictions of the 
time-dependent volume fraction transformed by 
isothermal annealing and computer simulations of 
the transformation of a finite ensemble of non-inter- 
acting spins on a lattice [19]. Those results for the 
volume fraction transformed are reproduced in Fig. 2 
for a scaled time, k, 

k = ut/R , (21 ) 

where u is the growth velocity of the cluster and R 
is the particle radius. For two nuclei per particle 



Fig. 2. The fraction transformed as a function of reduced time, k, 
for (a) one, (b) two and (c) six crystallites in a particle. The 
simulation results were obtained on a spherical lattice of radius 50 
units. The uncorrected JMAK (Eqs. (6) and (7)), finite-size cor- 
rected (FS) (Eq. (25) and large N (Eqs. (19) and (20)) solutions 
are shown. The Finite-size corrected solution provides the best fit 
to the lattice simulation. Taken from Ref. [19]. 
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(N = 2), the volume fraction transformed predicted 
from Eqs. (16), (17) and (20) fits the lattice simula- 
tion well to x = 0.5, but diverges markedly above 
that value, approaching a value less than complete 
transformation. As we will show, this incorrect 
asymptotic behavior is the origin of the predicted 
double peak in the DSC scan. The agreement be- 
tween the simulation and predicted results improves 
quickly with increasing number of nuclei per parti- 
cle, with the asymptotic value better approaching the 
correct value of 1. The simulation results and the 
predicted behavior from Eqs. (16), (17) and (20) are 
in almost perfect agreement for N > 6. Because of 
the low nucleation rates in silicate glasses, it is easily 
possible to produce powders of sufficiently small 
size that the number of growing clusters per particle 
is less than 6. To accurately model the transforma- 
tion behavior, then, it is desirable to extend the 
analytical expression for finite-size effects to those 
cases. 

For N spherical crystallites of the same size grow- 
ing in a spherical particle, the approximate expres- 
sion developed in Section 2.2.1 for the transformed 
volume fraction, taking account of overlap, can be 
rewritten as (i.e., combining Eqs. (14), ( 16), (17) and 
( 20 )) 

,v= 1 -exp[-yv(K 3 -^K 4 + ^/< h )]. (22) 

As shown, this provides an excellent description for 
x(f) when the number of clusters per particle, N, is 
greater than 5 to 10. To extend it to smaller sizes, it 
is useful to consider the extreme case of N — 1, for 
which Eq. (22) provides a poor description (Fig. 2). 
For the case of one crystallite growing in a spherical 
particle, there is no overlap, and Eqs. (14) and (16) 
may be used directly (recasting here for the scaled 
time), 

X = K } ~ j^K 4 + 3^K\ (23) 

which can be expanded to the same form as Eq. (22): 

- 1 -exp[-(K 3 -> J + S«‘-iK 7 +...)l. 

(24) 

Surprisingly, even for the extreme case of N — 
Eqs. (22) and (24) are identical up to the k 6 term. 
Based on this, Levine et al. [19] demonstrated that 
Eq. (24) can be modified simply to describe finite- 


size effects for all particle sizes, making the coeffi- 
cient of the k 6 term a function of the number of 
crystallites per particle and dropping higher-order 
terms: 


x = 1 — exp 



3 

K 



m( N) 
32 



(25) 


Fitting the predictions of Eq. (25) to the results of 
the lattice simulation gave 


m( N ) = max(5.26 - 0.26/V, 1.0). (26) 

The quality of the fits from Eq. (25) are shown in 
Fig. 2; the computed values for x(/) are virtually 
indistinguishable from the lattice simulation results. 
Two additional points should be made: (I) the x(t) 
computed from the JMAK solution becomes better 
as N increases, although it does so slowly (only 
approaching within 1% of the correct value at half- 
transformation when N = 42,000 [19]) and (2) Eq. 
(25) is valid even for a single nucleus growing in a 
particle, where overlapping volume need not be con- 
sidered. Eq. (25) then provides a simple expression 
for taking account of finite size transformation ef- 
fects over an astonishingly large range. As was 
demonstrated previously [19], Eq. (25) can also be 
modified easily to take particle shape into account by 
substituting for k an effective value, k c : 



(27) 


where V c is the extended volume of a growing 
crystallite and V is the volume of the particle. Eq. 
(27) was used to convert the transformation data 
computed for an infinite sample size (cf. Section 
3.1.). Since for lithium disilicate glass the nucleation 
and growth regions are w ? ell separated in tempera- 
ture, all crystallites will be approximately the same 
size. The reduced time, k c , is then directly obtained 
from the extended volume fraction transformed for 
an infinite system x*(r): 


<(<) 

N 


1/3 


(28) 


An analysis of the asymptotic value of Eq. (22) 
reveals the origin of the double peaks predicted in 
the DSC calculations when following the method 
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outlined in Section 2.2.1. That development is only 
valid for k < 2, corresponding to the maximum time 
required for a single cluster to transform the particle. 
From Eq. (22) that time corresponds to x = 1 — e lV , 
the asymptotic value observed in Fig. 2, showing 
why the agreement between the predicted x(t) and 
the lattice simulation results improves rapidly with 
increasing N. As k increases without bound, Eq. 
(22) will eventually approach 1, although this is 
outside the region of its validity. The first predicted 
peak, then, corresponds to the marked decrease in 
the rate of transformation corresponding to the first 
asymptote; the second peak arises because of the 
exponentially increasing growth velocity with tem- 
perature. carrying the expression to the second, al- 
though non-physical, asymptote. As expected, since 
Eqs. (25) and (26) give the correct asymptotic behav- 
ior (Fig. 2), only single DSC/DTA peaks are pre- 
dicted. 

2.2.3 . Surface crystallization 

As demonstrated in our companion paper, surface 
crystallization is appreciable in lithium disilicate 
glass. Owing to the increase of the surface-to-volume 
ratio, it can dominate the crystallization behavior of 
the small particles generally used in DSC/DTA 
experiments. The magnitude and temperature depen- 
dence of the surface nucleation rates are required to 
model this process, although for a polymorphic 
transformation, the crystal growth rate should be the 
same as that used for internal growth. Optical mi- 
croscopy showed that in the undoped or non- 
nucleated glasses, a thin, complete, surface layer 
grew toward the particle center, indicating that the 
surface nucleation rate is extremely high in this 
temperature range. For the purposes of modeling, 
therefore, we ignore the nucleation step and calculate 
only the rate of transformation due to a spherical 
growth front originating from the surface. 


Assuming a constant growth velocity, it. the vol- 
ume fraction transformed by surface growth alone is 


R - at 
R 


(29) 


Experimentally, surface and volume crystalliza- 
tion proceed simultaneously, making it necessary to 
account for overlap between the two processes as 

*(') =-*, + •*„(• --O- (30) 

Here jc v . is the fraction transformed from volume 
nucleation and growth. This approach is similar to 
previous treatments [8,9]. 


2.3. Input parameters 


Several parameters were required to model the 
experimental data for the crystallization of lithium 
disilicate glass. The free energy was calculated as a 
function of temperature, T [20]: 

AG — a^ + r/j T + a-, T~ T afT^ . ( 3 1 ) 

A temperature dependent interfacial energy was ob- 
tained by fitting the measured steady-state nucleation 
rates [21] 

<x = <x 0 + <t,7\ (32) 


The viscosity of the glass was assumed to follow a 
Fulcher- Vogel form: 


V = Vo ex P 


T-T 


(33) 


where tj 0 , and T a are estimated from experimental 
data [22]. The value for r| 0 was refined by fitting the 
measured transient times for homogeneous nucle- 
ation [21], keeping the temperature dependence fixed 
at the value obtained from the viscosity data. The 


Table I 


Gibbs free energy (J/mol) 

Ref. [20] 

a {) = 48045 

«, = -36.8! 

a 2 = 0.005607 

Interfacial energy (J/nr ) 

Ref. [21] 

c7 () = 0.094 

(r, = 7 x nr- 


Viscosity ( P) 

Ref. [22] 

7f kJ = 0.0363 

r*~: 

II 

T a = 460 

Atomic volume (10 -6 nr /mol) 

Ref. [23] 

r = 61.2 



Jump distance (A) 

Ref. [21] 

A = 4.6 
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diffusivity was calculated from the Stokes-Einstein 
relation 


D = 


kyj 

3itar} 


(34) 


where a is a characteristic distance of order of the 
atomic diameter. The values of the parameters used 
are listed in Table 1 . 

The value of C in Eq. (5) was determined by 
requiring agreement between the calculated and mea- 
sured [15,24] values for the macroscopic growth rate 
in lithium disilicate glass. The value, C ~ 4.8, so 
determined was held fixed for all calculations pre- 
sented here. 


3. Results and discussions 

The numerical approach was tested by comparing 
model predictions for DTA/DSC studies with exper- 
imental data from lithium disilicate (LS 2 ) glasses. 
Glasses prepared with different quench rates have 
different numbers of quenched-in nuclei and differ- 
ent cluster-size distributions. The simulated glasses 
used in these studies were quenched on the computer 
following a procedure described elsewhere [14]; un- 
less stated otherwise, a quench rate of l°C/s was 
assumed. Possible quench-rate effects on the atomic 
mobility are not considered. The simulated as- 
quenched glasses were annealed on the computer, 
under both isothermal and non-isothermal conditions, 
and scanned to higher temperatures at different rates 
to simulate a DSC/DTA scan. The results are pre- 
sented in this section. 

3.J. Particle size effects 

As discussed in our companion paper [1], the 
DSC/DTA peak parameters are a strong function of 
the particle size. To investigate the validity of the 
numerical model, DSC/DTA peaks for the as- 
quenched glasses were computed at a scan rate of 
15°C/min for several particle radii: 17.5 + 5 jam, 
45 + 8 jam, 120 + 30 fxm, 231 + 19 jam, and 510 + 
85 (am. These size ranges are the same as those used 
in the experimental measurements [1]. A uniform 
size distribution was assumed. 


3. }.}. Spherical particles 

Given an ensemble of small particles, if the nucle- 
ation rate is low it is likely that some particles will 
never develop internal nuclei. The statistical nature 
of this nucleation step was treated following the 
procedure discussed in Section 2.2. 

The results from calculations for spherical parti- 
cles undergoing finite-size volume (FS) and surface 
crystallization are compared with experimental data 
in Fig. 3, showing the peak temperature, T p (a) and 
the peak height (b) as a function of particle size. For 
comparison, the predicted behavior for the surface 
(Eq. (29)) and volume (FS, Eq. (25)) contributions 
and the asymptotic values for the infinite sample are 
also shown. There is little difference between the 
surface and surface + volume calculations of the peak 
temperature; both agree well with the experimental 
data. Although the magnitude of change in the peak 
temperature predicted from internal nucleation alone 



Fig. 3. Effect of particle size on (a) the DSC/DTA peak tempera- 
ture and (b) the relative peak height (height normalized to that of 
the 35 p.m particles). (Note that the upper part of the split vertical 
axis for (b) is logarithmic while the lower part is linear.) Calcu- 
lated results for an infinite sample, finite-size corrected (FS), 
surface crystallization and a combination of FS and surface cor- 
rections are shown. All calculations were made for spherical glass 
particles that w'ere quenched at l°C/s. 
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(FS) disagrees with the data, the trend with particle 
size is correct. Both the magnitude and direction of 
the change are incorrectly predicted for the peak 
height, however. This strong effect is a result of the 
low volume nucleation rate in these glasses; it re- 
flects the statistical absence of nuclei in a sizable 
fraction of the small particles. Without surface crys- 
tallization, those particles do not crystallize and hence 
do not contribute to the DSC/DTA signal. Larger 
particles have a higher probability of developing at 
least one nucleus. Consequently, a larger volume 
fraction will transform, resulting in a larger signal. 
This also explains why the peak temperature rises 
slightly above the infinite particle value for interme- 
diate particle sizes. The excellent agreement between 
the experimental data for small particles and the 
calculation based only on surface crystallization, and 
the near indistinguishability of the surface and sur- 
face + volume calculations for these small particles 
reflects the decreased importance of internal growth 
due to the large surface-to-volume ratio. Only at the 
largest particle sizes do the surface and surface + 
volume calculations begin to differ; both however 
follow the general trend of the data. 

These calculations show a broad minimum in the 
DSC/DTA peak height near particle diameters of 
400-600 fim, which becomes more prominent when 
the calculations are carried out to larger particle sizes 
(Fig. 4(a)). Three regions can be defined as a func- 
tion of particle size; they reflect behavior arising 
from a decreasing surface-to-volume ratio with in- 
creasing particle size. Region 1, occurring before the 
minimum in peak height, is dominated by surface 
crystallization. Both surface and internal crystalliza- 
tion are important in region 2 immediately following 
the minimum and extending to the point where the 
curve has nearly reached its asymptotic (infinite 
volume) value. Volume nucleation and growth is 
dominant for the large particle sizes of region 3. A 
similar, although less dramatic, crossover behavior is 
also predicted for the peak temperature (Fig. 4(b)). 
The continuing increase in peak temperature with 
particle size for the surface calculation simply re- 
flects the decreasing surface-to-volume ratio. The 
strong effect predicted for the DSC /DTA peak height 
makes it a potentially useful parameter for ascertain- 
ing at what size volume nucleation begins to domi- 
nate, allowing a rough assessment of the importance 



Fig. 4. Calculated results for finite-size corrected (FS) plus surface 
crystallization and an infinite sample (a) for the DSC/DTA 
relative peak height, / , (height normalized to that of the 35 pun 
particles) and (b) the peak temperature. Three distinct regions are 
evident in the behavior of the peak height: surface crystallization 
doininates in zone 1, surface and volume crystallization contribute 
in zone 2, volume crystallization dominates in zone 3. The glass 
was quenched at l°C/s. 

of surface crystallization. Although the predicted 
minimum in the peak height is not observed experi- 
mentally, the data were not taken for sufficiently 
large particle sizes that the effect becomes pro- 
nounced. 

3.1.2 . Ellipsoidal particles 

For simplicity, spherical particles were considered 
for most calculations. Optical and electron mi- 
croscopy studies, however, show that the particles 
have a range of shapes as well. To investigate possi- 
ble effects due to this, changes in the DSC/DTA 
peak profiles were computed for ellipsoids. Fig. 5 
compares the predictions for the ellipsoidal and 
spherical particles as a function of particle size. 
While the functional dependence of the peak temper- 
ature on particle size is virtually unchanged, the 
maximum in the peak width (not shown) and the 
minimum in the peak height (Fig. 5(b)) are shifted to 
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Fig. 5. The calculated effect of particle shapes on (a) the 
DSC/DTA peak temperature and (b) the relative peak height. / . 
The legend indicates the aspect ratios of the ellipsoidal particles 
considered. The glass was quenched at l°C/s. 


larger cluster sizes and are more shallow for the 
more eccentric particles. 

3.1.3. Heterogeneous nucieation 

Support for the validity of these computer predic- 
tions is found in the changes reported in our compan- 
ion paper [1] for the DSC/DTA peak parameters 
with the density of nucleating agents in heteroge- 
neously nucleating glasses. The peak height de- 
creased monotonically with increasing particle size 
for pure lithium disilicate glass; it initially decreased 
but subsequently increased dramatically with increas- 
ing particle size in glasses that had been heavily 
doped with Pt. Using thermodynamic and kinetic 
parameters determined from nucieation measure- 
ments of the Pt-induced heterogeneous nucieation 
rates [25], values for T p and peak height were com- 
puted as a function of spherical particle diameter for 
pure lithium disilicate glass and for glasses contain- 
ing 1 ppm and 5 ppm Pt impurities (Fig. 6(a) and 
(b)). The peak profile is dominated by surface crys- 
tallization for small particles and is independent of 


the amount of Pt added to the glass. For larger 
particles, however, a lower asymptotic temperature is 
obtained with increasing Pt dopant level, reflecting 
the accelerated transformation due to the larger num- 
ber of internal nuclei. Similarly, the peak height goes 
through a minimum that occurs earlier and is sharper 
in the more heavily doped glasses. Although the 
trends are predicted correctly, the calculated results 
disagree in detail with the experimental data; in 
particular the sharp increase experimentally observed 
in the heavily doped (5 ppm) glass is not reproduced 
and the predicted rise in peak height for the 1 ppm 
doped glass is not observed. 

The predicted doping dependence of the peak 
temperatures for ellipsoidal particles has the same 
form with increased particle size as that shown in 
Fig. 6(a), although the magnitudes are different. As 
expected, the effect on the peak height is more 
dramatic (Fig. 6(c)). Significantly better agreement is 
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Particle diameter (pm) 


Fig. 6. The predicted effect of Pt doping on (a) the peak tempera- 
ture and (b) the relative peak height. / for spherical particles. 
Figure (c) shows the predicted change in the relative peak height 
for ellipsoidal particles of aspect ratio (I, f, 1 ). The legend 
indicates the amount of Pt introduced by weight: the points 
indicate the experimental data: the lines represent the computed 
result. All calculations were performed for glasses quenched at 
rc/s. 
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obtained between the computer predictions for the 
ellipsoidal particles and the experimental data on the 
undoped and 1 ppm doped glasses. The heavily 
doped glasses, however, still show' a sharper increase 
than is predicted from the calculations. This suggests 
that the assumption of ellipsoidal particles alone still 
may not adequately capture the range of particle 
shapes. More detailed comparisons would require 
tighter control of the particle shape distribution of 
the glasses. 

3. 1.4. Effects of water 

Water is known to have a profound effect on the 
stability of silicate glasses [26], presumably due to an 
increase in the atomic mobility. This was evidenced 
by a difference in the experimentally observed 
DSC/DTA peak profile parameters for "dry’ glasses 
versus k wef glasses [1]. A lower transition tempera- 
ture and a higher peak height were reported in 
glasses that were kept in a humid environment for 
extended times. Since the diffusion of water into the 
sample occurs from the surface, the enhanced surface 
crystallization is expected to have the greatest impact 
on the transformation kinetics. The surface growth 
will only accelerate, however, to a depth correspond- 
ing to the diffusion distance of the water. 

To qualitatively model these data, a standard dif- 
fusion profile for one-dimensional diffusion of water 
into the particle was calculated for the time that the 
sample was placed in the most atmosphere (80 h at 
room temperature). No account was taken of water 
re-distribution during the subsequent DSC/DTA 
scan. Lacking data for the diffusion coefficient of 
water in lithium disilicate glass, measured values for 
water diffusion in fused silica [26] were used to 
estimate a room temperature diffusion coefficient, 
giving D ~~ 2.8 X 10 10 cnr/s. The nucleation rate 
in hydrated glasses is known to be higher by almost 
an order of magnitude over that in non-hydrated ones 
[27]. Assuming that this is entirely due to a change in 
the mobility, a multiplicative constant for the growth 
velocity was taken to scale linearly with the water 
concentration, taking on values between one for the 
dry glass and ten for the fully hydrated glass (at the 
surface of the particle). Fig. 7(a) and (b) show' the 
peak temperature and the peak height computed as a 
function of particle size. They are in fair agreement 
with the experimental data, showing a lower transi- 



Particle diameter (gm) 

Fig. 7. The predicted effect of water on (a) the peak temperature 
and (b) the relative peak height, / p . The points represent the 
experimental data; the lines represent the calculated results. Water 
treatment: dry samples (solid line, ■ ); surface hydrated sample 
where water diffused in from the surface (dotted line. •); fully 
hydrated sample where water was assumed uniform throughout 
the volume (dashed line). All calculations were performed for 
spherical glass particles quenched at l°C/s. The experimental 
error is of order the si/e of the data points. 

tion temperature and a higher peak height in the 
hydrated glasses for the smaller particles. For the 
smaller particles, the calculations predict a difference 
in the peak temperatures and peak heights for the 
hydrated and dry glasses that is not observed experi- 
mentally. The reasons for this are not known; as for 
the Pt doped glasses, it may indicate a distribution of 
particle shapes different from those modeled here. 
Because the thickness of the hydrated layer ex- 
pressed as a fraction of the particle diameter de- 
creases with increasing particle size, the computed 
DSC/DTA parameters for the wet and dry glasses 
approach one another, and are in better agreement 
with the data. For comparison, the calculated 
DSC/DTA parameters for glasses that are fully 
hydrated throughout the volume are also shown in 
Fig. 7. Since both the surface and the volume growth 
velocities and the volume nucleation rate are in- 
creased by a factor of ten for all particle sizes, the 
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parameters remain different for all particle sizes. 
While these calculations suggest that DSC/DTA 
data could be used to make rough estimates of the 
diffusion coefficient of water in glasses, a quantita- 
tive comparison would require an iterative calcula- 
tion with a varying diffusion coefficient, modeling 
the change in both growth velocity and the nucle- 
ation rate as a function of distance from the surface. 

3.2. N on-isothermal nucleation treatments 

We have demonstrated previously that DSC/DTA 
scans following preannealing treatments can be used 
to estimate the temperature range for significant 
nucleation [4,7]. Here we compare model predictions 
for these studies with experimental data. Simulated 



<f> n (C/min) 

Fig. 8. The calculated effect on the DSC/DTA peak parameters 
of scanning through the nucleation /.one at different rates: (a) the 
peak temperature (h) the peak height, / p , and (c) the peak 
half-width. Results are shown for infinite samples (dot-dot-dash 
line), finite-size corrected (PS) (solid line), surface crystallization 
(dashed line) and a combination of FS and surface corrections 
(dotted line). Experimental data are also shown: (#) DSC data 
and ( □ ) DTA data. All calculations were made for glass particles 
uniformly distributed between 425-500 |xm in diameter that had 
been quenched at 0. PC/s. The experimental error is of order the 
size of the data points. 
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as-quenched glasses were scanned first through the 
nucleation zone, from 400°C to 500°C. at rates (0 N ) 
between 2 and 15°C/min and were subsequently 
scanned through the crystallization peak. The crystal- 
lization scan rate. ( P C . was typically 15°C/min. 

Fig. 8 compares the predictions of the computer 
model for the key DSC/DTA peak parameters with 
experimental data as a function of the scan rate 
through the nucleation zone. Spherical particles of 
simulated glasses quenched at 0. l°C/s w ith an aver- 
age diameter of 462 jxm were assumed. The results 
from two independent DSC/DTA measurements 
made in different laboratories on these glasses are 
provided. The scatter between these data provides 
some indication of the expected experimental error. 
They indicate that the peak temperature and the peak 
height are the more robust parameters, better suited 
for quantitative comparison with model predictions. 
For comparison, calculated values for an infinite 
sample, for finite-size corrections only, for surface 
crystallization alone, and for combined surface and 
volume crystallization for finite-size particles are 
shown. In all cases, reasonable agreement with ex- 
perimental observation is found. As expected, T p and 
A T p increase with increasing scan rate through the 
nucleation zone while / p decreases. As the scan rate 
through the nucleation zone is decreased, the time 
for nucleation increases, resulting in a greater popu- 
lation of nuclei, causing an increase in the transfor- 
mation rate. Although the agreement between the 
predicted and measured peak temperature is less 
good, it is still reasonable. It should be emphasized 
that no fitting parameters w p ere adjusted for these 
calculations; the input data are those obtained from 
experimental measurements of the time-dependent 
nucleation rate and the macroscopic growth velocity 
as a function of temperature. 

The calculated peak parameters for two particle 
sizes, 462 |xm and 1020 (xm, of a glass quenched at 
0. l°C/s are compared with experimental data in Fig. 
9. Spherical and ellipsoidal particles were consid- 
ered. In all cases, the calculated results predict the 
correct trends, although they disagree quantitatively 
with the measurements. For the spherical particles 
(Fig. 9(a), (b) and (c)), the computed values for T p 
show a larger change than is observed in the data, 
while the computed values for 7 p and A T p are less 
than is observed. The agreement with experimental 
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Fig. 0. The effect of scanning through the nucleation zone at different rates on the DSC/DTA peak parameters. Experimental data from 
DSC (■) and from DTA (□) for particles with diameter 425-500 pm; DSC data from particles with diameter 850-1170 pm (•). 
Calculations for spherical particles with diameter 425-500 pm (solid line) and 850-1 170 pm (dashed line): (a) the peak temperature (b) the 
peak height, / p , and (c) the peak half-width. Corresponding results for ellipsoidal particles of eccentricity (1. y, y) are shown in (d). (e). and 
(f) respectively. All glasses were quenched at 0. l°C/s. Calculations are shown for surface plus finite-size corrected volume crystallization. 
The experimental error is of order the size of the data points. 
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Fig. 10. The effect of scanning through the nucleation zone at different rates on the DSC/DTA peak parameters. Experimental data from 
DSC (■) and from DTA (□) for particles with diameter 425-500 pm; DSC data from particles with diameter 850-1170 pm (•). 
Calculations for spherical particles with diameter 425-500 pm (solid line) and 850-1 170 pm (dashed line): (a) the peak temperature (b) the 
peak height. / p . and (c) the peak half-width. Corresponding results for ellipsoidal particles of eccentricity ( 1, y, y) are shown in (d), (e). and 
(f) respectively. All glasses were quenched at l c C/s. Calculations are shown for surface plus finite-size corrected volume crystallization. 
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Fig. 1 1. Calculated relative peak heights. / p . lor glasses that had 
been scanned through the nuc leal ion zone (solid line) and for 
glasses held isolhermally for a corresponding amount of time at 
the peak temperature of the steady-state nuclealion rate (dashed 
line) compared with corresponding experimental data ((• ) scanned 
glasses and (O) isothermallv treated glasses). Calculations in (a) 
were made assuming time-dependent nucleation: steady-state nu- 
cleation was assumed in (b). All glasses were quenched at 0. 1 C /s. 
The experimental error is indicated by the error bars. 


data is significantly improved by assuming ellip- 
soidal particles instead of spherical ones (Fig. 9(d), 
(e) and (0). although the predicted change in peak 
height as a function of particle size is smaller than is 
experimentally observed in both cases. 

Since the experimental quench rates for the glasses 
were unknown, these calculations were repeated on 
simulated glasses that were quenched from the melt 
at rates of l°C/s and 10°C/s. The results for the 
10°C/s quenched glasses were almost indistinguish- 
able from those for the l°C/s glasses, which are 
shown in Fig. 10. The peak temperatures are fit 
slightly better by the glasses quenched at 0. l°C/s, 
although the differences with calculations for glasses 
quenched at l°C/s are small. 

Fig. 11 shows the computed DSC/DTA peak 
heights for glasses that were either scanned at differ- 
ent rates through the nucleation zone (400 o C-500°C). 


or were annealed isothermal ly (following the meth- 
ods in Refs. [6,7]) at the nucleation peak temperature 
for a length of time equal to the total time required 
to scan through the range of significant nucleation 
(425°C-500°C). As observed experimentally, the 
calculated peak heights for the two annealing treat- 
ments are identical, whether time-dependent (Fig. 
11(a)) or steady-state nucleation (Fig. 11(b)) rates 
are assumed, in agreement with the experimental 
results. The reasons for this unexpected behavior are 
not obvious and are discussed further in the follow- 
ing sections. The experimentally observed changes in 
peak height with increasing scan rate and decreasing 
annealing time are sharper than predicted from the 
steady-state calculations, indicating a stronger depen- 
dence on the production of nuclei than is expected 
from the temperature dependence of the steady-state 
nucleation rate alone. Importantly, the experimental 
and calculated results are almost indistinguishable 
when time-dependent nucleation effects are included. 

We have shown previously that the cluster distri- 
bution and hence the nucleation rate of as-quenched 
glasses is considerably decreased from the steady- 
state value [6,14]. Computer calculations of the sub- 
sequent nucleation behavior on non-isothermal heat- 
ing predicted that the nuclealion rate would rise 
above the steady-state value, peaking at a tempera- 
ture higher than the peak of the steady- state rate [6]. 
Similar behavior is evident following the multi-step 
annealing treatments considered here. Fig. 12(a) 
show's the calculated nucleation rate at a cluster size. 
n = 310 (well above the largest critical size, n , in 
the temperature range over which the transformation 
occurs) as a function of temperature for glasses 
scanned through the nucleation zone at <i\ = 4 and 
0.5°C/min and for glasses annealed isothermallv for 
corresponding times of 18.73 and 1 50 min. In all 
cases, the nucleation rate rises above the steady-stale 
rate at a temperature that is higher than that of the 
peak in the steady-state rale. Interestingly, two peaks 
in the nucleation rate are observed when 0 n = 
0.5°C/min. The origin of this complex behavior is 
found by an examination of the temperature depen- 
dence of the cluster distribution. From Eq. (3) for the 
large clusters considered, where 
N n M f , the nucleation rate is l n! = N„ ,( k~ — k n ). 
which is in turn approximately A r nI k f \ since the 
forward rate is much greater than the backward rate 
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Fig. 12, (a) The calculated nucleation rate. /, as a function of 
temperature for glasses annealed in the nucleation /one: (b) the 
corresponding total number of nuclei. A'. The solid and dashed 
lines are for samples scanned through the nucleation zone at 
</>, = 4 and ().5"C/min. respectively: the dotted and dot-dot-dash 
lines are for samples annealed at the peak nucleation temperature 
for 18.755 and 150 min respectively. For comparison, the steady- 
state nucleation rate is also given. The nucleation /one (400- 
500°C) and a portion of the crystallization /one (500°C through 
the crystallization peak) are indicated. 

there. The computer calculations, of course, follows 
Eq. (3) exactly; these simplifications are only made 
here to facilitate the discussion. During a DSC/DTA 
scan to higher temperature, therefore, the observed 
nucleation rate will result from a competition be- 
tween the density. N n r which is evolving toward a 
temperature dependent steady state distribution, and 
the net rate, which is becoming larger with increas- 
ing temperature because of the increasing atomic 
mobility. Considering first the case where <P n = 
0.5°C/ min, reaches its maximum value near 

480 c C (Fig. 13), explaining the first peak in the 
nucleation rate. For this slow scan rate, the time 
spent at temperatures within the nucleation zone is 
sufficiently long that considerable relaxation toward 
the steady-state density occurs. Since tV, l(U(0 de- 
creases for all higher temperatures, the second, 
smaller, peak arises from the competition between 


the abrupt change in the rate of relaxation (and hence 
the rate of decrease of /V M(}7(f) ) when the scan rate 
is changed from <£> n to 0, , and the increased rate of 
change in k t \ resulting from the more rapid increase 
in temperature for ( P C . The effect becomes less dra- 
matic as 0 n approaches </>, finally disappearing 
when 0 n = 0. For 0 n = 4°C/min, only a kink is 
evident at the temperature for which the rate is 
increased to 0 C = 15°C/min, reflecting the sudden 
change in rate constants. Now, the peak in the 
nucleation rate corresponds closely with that of the 
maximum in cluster density, pushed to higher tem- 
peratures because of the less complete relaxation 
occurring within the nucleation zone. This difference 
in degree of relaxation of the cluster density is also 
evident in calculated behavior for glasses annealed at 
the peak nucleation temperature (450°C) for different 
times and subsequently scanned to higher tempera- 
tures at l5°C/min. The jump in the nucleation rate 
near 450°C for the glass annealed for 150 min is a 
result of the higher density. The more gradual change 
for the glass annealed for less time reflects the 
smaller degree of relaxation from the as-quenched 
distribution. At temperatures below 450°C, the nu- 
cleation rate is independent of the annealing treat- 
ment. The observed temperature dependence reflects 
that of k' u ; the cluster density remains virtually 
unchanged from the as-quenched value because of 
the low’ rate of relaxation. 



Fig. 13. 13A portion of the cluster density for different tempera- 
tures for tp n = ().5°C/min and <I{ = 15°C/min. The cluster si/.e at 
which the nucleation rate is measured (Fig. 12(a)). h = 310, is 
indicated. 
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Fig. 12(b) shows the corresponding number of 
nuclei produced as a function of temperature. Glasses 
annealed more slowly (or for longer time) in the 
nucleation zone contain more nuclei, although the 
reason for this is far more complicated than could be 
ascertained from steady-state considerations. For ex- 
ample, for the non-isothermal case where 0 n = 
0.5°C/min, most of the nucleation occurs in the 
nucleation zone. Approximately the same number of 
nuclei are produced in the corresponding case of a 
150 min isothermal anneal at the peak temperature of 
the steady-state nucleation rate, although they are 
nucleated at temperatures above the upper limit of 
the nucleation zone. Similar considerations apply for 
the higher values of <f> n and the corresponding 
isothermal treatments. Anneals in the nucleation 
zone, then, not only result in the production of 
nuclei, but also cause differing degrees of relaxation 
of the cluster distribution toward the steady-state 
distribution. Although there may be no significant 
production of nuclei at temperatures near the 
steady-state peak temperature, those clusters that are 
‘swept through’ the distribution result in a larger 
nucleation rate at higher temperatures. It is only 
fortuitous that in this glass the isothermal anneals 
for three-quarters of the time required to scan 
through the nucleation zone produce a similar num- 
ber of nuclei in the scanned samples . 

These results further demonstrate the role of 
time-dependent nucleation phenomena in determin- 
ing the stability and transformation behavior of 
glasses and other phases. They also show that simple 
interpretations of DSC/DTA non-isothermal trans- 
formation data are not possible. The success of the 
numerical model, however, indicates that it can be 
used to model these data to obtain quantitative ki- 
netic information. 

3.3 . Analysis of peak profiles 

The differences in peak parameters discussed in 
the previous sections indicate that it should be possi- 
ble to distinguish between various crystallization 
mechanisms based on a qualitative inspection of the 
peak shapes. For illustration, the computed DSC 
peak profiles for a heating rate of 15°C/min of a 
sample consisting of 50 (xm diameter particles are 
shown in Fig. 14. These glasses were quenched from 



Fig. 14. The computed DSC/DTA peaks for 50 jjim diameter 
glass particles: (a), (d) infinite particle size; (b), (e) spherical 
particles; (c), (f) ellipsoidal particles with eccentricity ( 1 , y. y). 
Figures (a), (b) and (c) were for as-quenched glasses; (d), (e) and 
(f) were for glasses that were scanned at l°C/min between 
400-500°C followed by a 15°C/min scan from 500-900°C. All 
calculations were performed for glasses quenched al l c C/s. 


the melt at l°C/s and had no preannealing treat- 
ments prior to the DSC/DTA experiment. The peak 
shape for the infinite sample (Fig. 14(a)) is qualita- 
tively different from that obtained when finite parti- 
cle-size effects and surface crystallization are in- 
cluded. The transformations for finite-sized spherical 
(Fig. 14(b)) and ellipsoidal (Fig. 14(c)) particles 
occur at lower temperatures and are spread out over 
a larger temperature range. Such differences could be 
interpreted by an appropriate model of the growth 
velocity, however. What might be interpreted experi- 
mentally is the degree of peak asymmetry. Peaks for 
the small particles show a much longer tail on the 
low temperature side of the peak. A comparison with 
the peaks from surface and volume crystallization 
alone shows that the peak shape for such small 
particles is dictated largely by surface crystallization. 
The predicted peak for ellipsoidal particles looks 
very similar to that for spherical particles, although 
the drop-off at high temperatures is sharper, reflect- 
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l ig. 15. The computed DSC/DTA peaks for 500 p.m diameter 
glass particles: (a), (d) infinite particle si/e: (b). (c) spherical 
particles; (c). (D ellipsoidal particles with eccentricity (I. , ). 

1 igurcs (a).(h) and <c) were for as-quenched glasses; (d). (e) and 
(f) were for glasses that were scanned at I’C/min between 
400-500 (' followed by 15 C /min scan from 500-000 C'. All 
calculations were performed for glasses quenched at I'C/s- 

ing the greater importance of surface crystallization 
due to the shorter distance between surfaces in two 
of the three dimensions. Fig. 1 4( d )— (f ) show the 
predicted peaks for the infinite, and finite-sized 
spherical and ellipsoidal particles after the glass has 
been heated through the nueleation zone at l°C/min. 
That the peaks predicted for the smaller particles are 
virtually unchanged by the preannealing treatment 
reflects the dominance of surface crystallization, 
which is unaffected by annealing in the nueleation 
zone, since the surface nueleation step is ignored 
here. The peaks for the infinitely large particles are 
shifted to lower temperatures because of the in- 
creased number of nuclei resulting from the prean- 
neal. 

Fig. 15 shows the predicted peaks for larger, 500 
(urn diameter, particles of as-quenched glasses and ol 
glasses that have been heated at TC/min through 
the nueleation zone. Taking both surface and volume 
effects, the predicted DSC/DTA peaks for particles 


of the as-quenched glasses are still more asymmetric 
than for an infinite glass, showing a longer tail on 
the low temperature side of the peak. They are 
displaced to higher temperatures than for the smaller 
particles (Fig. 14(b) and (c)), reflecting the decreased 
importance of surface relative to volume nueleation. 
Annealing within the nueleation zone changes the 
DSC/DTA peak significantly for both the spherical 
and ellipsoidal particles. The peak temperatures are 
shifted down and the peak heights are increased, 
becoming more similar to the behavior of the infi- 
nite-sized particles. The increasing importance of 
volume nueleation and growth is also evidenced by 
the more symmetric peak shapes. 

4. Conclusions 

In conclusion we have presented a model for 
simulating polymorphic crystallization under non- 
isothermal and isothermal annealing conditions for 
particles of finite size. Surface and volume nuele- 
alion and growth, time-dependent nueleation rates 
and cluster-size dependent growth velocities are in- 
cluded. Predicted DSC/DTA scans of lithium disili- 
cate glasses that had been heated at different rates 
although the temperature range where the steady-state 
nueleation rate is significant w r ere compared with 
experimental results. The calculations were made 
using independently measured values for the nucle- 
alion and growth rates, the temperature-dependent 
free energy difference between the glass and the 
crystal phases and the viscosity of the glass. Calcula- 
tions were made as a function of particle size and 
particle shape and compared with the experimental 
data. Glasses quenched at different rales were also 
investigated. 

Calculated peak profiles were in very good agree- 
ment with measured data, given that no parameters 
were adjusted from their known values to improve 
the fits. The DSC/DTA peak height was shown to 
he sensitive to the contributions of the surface and 
volume and to the particle shape. A minimum in the 
peak height was predicted near particle radii of 200 
pan, due to the increasing importance of volume 
over surface crystallization with increasing particle 
size. Evidence for this minimum was found by exam- 
ining the crystallization behavior of glasses as a 
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function of Pt dopant level, a known heterogeneous 
nucleating agent that increases the contribution of 
internal growth. Future experiments on glass parti- 
cles of carefully controlled size and shape will pro- 
vide a more quantitative test of the model* 

The good agreement strongly supports the model. 
It can be used to more accurately model DSC/DTA 
data than has been possible previously. Further, the 
numerical model provides a more complete insight 
into the fundamental processes of the devitrification, 
enabling better control of processing parameters and 
allowing estimates to be made of key kinetic parame- 
ters. The effects of time-dependent nucleation on 
phase stability and the transformation behavior and 
the ability of our numerical technique to model these 
quantitatively was demonstrated. The numerical 
model can also be used to design new experimental 
methods for analyzing calorimetric data. For exam- 
ple, these calculations demonstrate that if surface 
crystallization is known to be important in a system, 
the size below' which it is the dominant mechanism 
can be determined from a plot of the height of the 
DSC/DTA peak versus particle size for particles of 
knowm shape. In those cases, the growth rate can be 
estimated as a function of temperature from one 
DSC/DTA experiment by assuming surface crystal- 
lization to fit the DSC/DTA peak. 
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